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Abstract. - Eftects of a coupling between the mechanical vibrations of a quantum dot placed 
between the two leads of a single electron transistor and coherent tunneling of electrons through 
a single level in the dot has been studied. We have found that for bias voltages exceeding a 
certain critical value a dynamical instability occurs and mechanical vibrations of the dot develop 
into a stable limit cycle. The current- voltage characteristics for such a transistor were calculated 
and they seem to be in a reasonably good agreement with recent experimental results for the 
single Ceo-molecule transistor by Park et al. ( Nature 407, (2000) 57). 



Introduction. - Nanoelectromechanics |0, D is a new, quickly developing field in con- 
densed matter physics. A coupling between strongly pronounced mesoscopic features of the 
electronic degrees of freedom (such as quantum coherence and quantum correlations) and de- 
grees of freedom connected to deformations of the material produces strong electromechanical 
effects on the nanometer scale. The mesoscopic force oscillations in nanowires [^[-^ observed a 
few years ago is an example of such a phenomenon. Investigations of artificiaUy-made nanome- 
chanical devices, where the interplay between single-electron tunneling and a local mechanical 
degree of freedom significantly controls the electronic transport, is another line of nanoelec- 
tromechanics P-[l5[|. For one of the nanomechanical systems of this kind, the self-assembled 
single-electron transistor, a new electromechanical phenomena - the shuttle instability and 
a new so-called shuttle mechanism of the charge transport were recently predicted in p^ . 
It was shown that a small metallic grain attached to two metallic electrodes by elastically 
deformablc links breaks into oscillations if a large enough bias voltage is applied between the 
leads. For the model system studied in it was also shown that a finite friction is required 
for the oscillation amplitude to saturate and for a stable regime of oscillations to develop. 

An essential assumption made in [ p^ is that the relaxation mechanisms present are strong 
enough to keep the electron systems in each of the conducting parts of the transistor in local 
equilibrium (as assumed in the standard theory of Coulomb blockade 1 16, ^). Such relaxation. 



which destroys any phase coherence between electron tunneling events, allows a description of 
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Fig. 1 - Model system consisting of a movable quantum dot placed between two leads. An efltective 
elastic force acting on the dot from the leads is described by the parabolic potential. Only one single 
electron state is available in the dot and the non-interacting electrons in each have a constant density 
of states. 



the whole electronic kinetics by means of a master equation for the occupation probabilities of 
the dot. It is clear that such an incoherent approach must fail if the size of the grain is small 
enough. Firstly, a decrease of the grain size results in a decrease of its mass and consequently 
in an increase of the frequency of the mechanical vibrations, which for nanometer-size clusters 
is of the order of 0.1 ITHz and might exceed the electron relaxation rate in the grain. 
Secondly, the electron energy level spacing in a nanometer-size grain S ^ £f/N, (TV is the 
number of electrons) can be of the order of 10 K and exceed the operational temperatures. In 
this situation the discreteness of the energy spectrum can be very important. What is more, 
a possibly strong tunneling-induced coupling between electronic states in the grain and in the 
leads, which results in large quantum fluctuations of the charge in the grain, must be included. 

For all of the above reasons a new approach is needed for a description of the electron 
transport through a nanometer-size movable cluster or quantum dot. The non-trivial question 
then arises whether or not the coherent electron tunneling through a movable dot causes any 
electromechanical instability. Such a question is of notable practical significance in view of the 
recent experiment by Park et al. , where the current through so-called single-Cgo-molecule 
transistors was measured and anomalies — including a few equidistant step-like features — in 
the I-V characteristics observed. This current steps were interpreted in as a manifestation of 
the coupling between electron tunneling and the center-of-mass vibrational degree of freedom 
of the molecule. 

Theoretical model. - To investigate the influence of the above electromechanical coupling 
we consider a model system consisting of a movable quantum dot placed between two bulk 
leads. An effective elastic force acting on the dot due to interaction with the leads is de- 
scribed by the parabolic potential presented in flg. |l|. We assume that only one single electron 
state is available in the dot and that the electrons in each lead are non-interacting with a 
constant density of states. At the same time we treat the motion of the grain classically. The 
Hamiltonian for the electronic part of the system is 

H = ^(e^fe - Ha)al^^.aak + ^d{t)c^ c + ^ Ta{t){a}^^c + c^aak) ■ (1) 

a,k a,k 

Here Tl^r = ri^fl exp{=Fx(t)/A} is the position-dependent tunneling matrix element, ed{t) = 
Co — ^x{t) is the energy level in the dot shifted due to the voltage induced electric field 
£/g = x^i X is a parameter characterizing the strength of the electrical field as a function 
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of the bias voltage V, e < is the electron charge, x measures the displacement of the dot, 
^ak creates an electron with momentum k in the corresponding lead, a = L, R is the lead 
index, creates an electron in the dot and A is the characteristic tunneling length jl2| . The 
first term in the Hamiltonian describes the electrons in the leads, the second — the movable 
quantum dot and the last term — tunneling between the leads and the dot. 

The evolution of the electronic subsystem is determined by the Liouville-von Neumann 
equation for the statistical operator /5(i), 

tdtm^[H,m]-, (2) 

while the center of mass motion of the dot is governed by Newton's equation, 

i + wlx = F/M . (3) 



Here wq = y^k/AI, M is the mass of the grain, fc is a constant characterizing the strength of 
the harmonic potential, F = — < dH/dx > and < • >= Tr{p{t)»}. The force F on the RHS 
of eq. is due to the coupling to the electronic subsystem and consists of two terms, 

F(i) = -zfG<(t,i) + 2A-i^(-irT„(t)Im[G<,(i,t)] , (4) 

where G^{t,t') = i{c^ {t')c{t)), G'^^,{t,t') = i{a'^^^.{t')c{t)) are the lesser Green functions, and 
a = (1) for the left (right) lead. The first term describes the electric force that acts on the 
charge in the dot; the second term is an exchange force, which appears due to the position 
dependence of the tunneling matrix elements Tl h. The force F depends on the correlation 
functions G^(t, t) and G'^i^{t, t), which can be computed exactly in the wide-band limit {pa = 
const) by using the Keldysh formalism fl^ j. 

Following the standard analysis jl^ we express the correlation function G'^^,{t,t) in terms 
of the Green functions of the dot, 

G<fe(i,t')= / At,TM{G^tM)9Uh,t') + G<{tM)glu{ti.t')] ■ (5) 

Here G^' is the retarded Green function of the dot and g^^^.^'' is the advanced (lesser) Green func- 
tion of the leads for the uncoupled system. The Dyson equation for the retarded (advanced) 
Green function has the following form: = ^''^'^^ + gri<^)j:ria)Qr(a) ^ ^j^g^g S''('^)(ii, ta) = 

J2c,k'^<='i*i)9alhti,t2)Ta{t2). In the wide-band limit E''^(")(ti, ts) oc (5(^1-^2) and this Dyson 
equation can be solved exactly. The lesser Green function is given by the Keldysh equation 
G< =G'■S<G^ where S<(ii, t2) = Eo,fe ^2)T„(t2). 

As a result we obtain a general expression for the force F of the form 

F{t)^Y.P-^j de/«(e)|£: |B„(e,t)|V2^^r„(0Re[B„(e,t)]| , (6) 

where ^ ^ 

B„{e,t) = -iJ diiT„(ti)exp|z^ dt2 e - £^(^2) + |, 

r{t) = 27r^^ paT^(t), Pa is the density of states in the corresponding lead and fai^) = 
[exp(/3(e — /ict)) -t- 1]^^. It is worth mentioning that eq. (||) is valid for arbitrary values of the 
tunneling matrix elements. 

An important question is now whether or not the mechanical stability of the transistor 
configuration is affected by coherently tunneling electrons. 
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Instability. - In order to investigate the stability of the equihbrium position of the 
dot we expand the RHS of eq. (||) to first order with respect to the displacement: F{t) = 
At'Dit — t')x{t'). After a Fourier transformation of the obtained equation we get the 
following dispersion equation for the frequency: 



^2 - {wl~wl) = -D^/M, 



(7) 



where 



A, 



Ra — 



de x ^ 

a 

(-1) 



A 



^afa{RaG^ + -w} ' 

(-1)" 



£\Go\ 



A 



-eo±«r/2]- 



^1 = Ea r„ / deURc[G+]/{7:XH'I), = 2^p„r2, T = To, G±(e) = [e^ 
and g = (Fl - rfl)/A. 

The criterion for instability is that the frequency u has a positive imaginary part. When 
T/iMw^X"^) < 1 and £:A/( ^w^o + T^) < 1 (the case of weak electromechanical coupling) 
we obtain that Im[u'] « —lm[Dwg]/{M'Wo) ■ From now on we consider only the case of weak 
electromechanical coupling. An analytical analysis can not be performed in the general case, 
but for a large symmetrically applied bias voltage and zero temperature one can show that 
Im[_Du,,-|] is negative: 

T rn 1^ ^^'^0 ^^^^ cc^ 
A^Wq+F-^) F 

The instability that follows from eq. (||) is in contrast with the behavior at low voltages 
where it can be shown that Im[Z3u,p] is positive. Therefore, a finite threshold voltage for 
the instability exists in the system. A simple expression for the threshold voltage can be 
found for a symmetric junction in the case of weak tunneling under the above conditions: 
eVc = 2(eo + wq). 

Under the condition of weak electromechanical coupling, the displacement x{t) of the dot 
can be represented in a harmonic oscillation form A{t)cos{wt + 4>{t)) with an amplitude A{t) 
and a phase (j){t) that slowly vary on the scale of w^^ . By averaging over the fast harmonic 
oscillations one can get the following equation for the evolution of the amplitude (see for 
example |p2|): 

dA^ _W{A^) 
dt MnwQ 

Here PF(A^) = <f dxF is the work done by the force F on the dot during one period of 
oscillation with a constant amplitude A. We conclude from eq. (^ that the stable regime of 
oscillations corresponds to W = and dW{A^) /dA^ < 0. Typical H^(A^-curves are depicted 
in fig H for the parameters taken from the experiment described in [Q. When an applied 
voltage is lower than the threshold voltage the work W is negative for all amplitudes. This 
implies that the equilibrium position of the grain is stable. When the voltage is higher than 
the threshold voltage the function W{A'^) is positive for < ^ < Ac{V) (which corresponds 
to a slow increase of the oscillation amplitude), negative for A > Ac{V) (which corresponds 
to a slow decrease of the amplitude) and equal to zero at A = Ac{V). This means that 
when the applied voltage exceeds the threshold value the amplitude of the oscillation slowly 
increases until it develops into a stable limit cycle with amplitude Ac{V). As one can see from 
fig. |[ the amplitude of the limit cycle is of the order of A ~ 0.1 nm and is considerably larger 
than the amplitude of zero-point oscillations of the grain, which for the Ggo molecule in the 
experiment & is approximately 3 pm. 
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Fig. 2 - The energy W pumped into the mechanical degree of freedom during one period of osciUation 
for different values of the applied voltage V: wq = 5meV, T — 0.13 meV, eo — 6meV, F = 1 /leV and 
Tr/Tl = 9. 



Current. - In the regime where the stable limit cycle oscillations have developed the 
vibration-induced inelastic tunneling of electrons can significantly contribute to the current. 
The time-averaged current through the system in the stable regime where oscillations of the 
limit cycle amplitude Ac(y) and frequency wq — ^JkjM have developed has the following 
form: 



2-kT 



dt de < 



Ba 



2e ^/Llm 



Bl 



(10) 



where x{t) — Accos{wot) and T = 2tt/wq is the period of oscillations. This expression is valid 
under the same general conditions as eq. (||). 

Under the condition that the effective level broadening T wo (where T — TJo{i2Ac/ X)) 
eq. (^^ can be simplified to 



eT r 



2m 
L 



(11) 



where /„,m = /a(eo + rnwo), r] = ^(f/wo)^ - A = -[{£/wo) + {-IY'/X]/t] and Jm are 

Bessel functions of the first kind. The typical TV curves are shown in fig. |3|, where we assume 
that voltage has been applied only to the left lead (according to the experiment 0). The 
main characteristic feature of all obtained TV curves is that they show only a few equidistant 
steps which are followed by step-less behavior of the curves. The distance between the steps is 
given by the vibrational frequency wq and their heights can vary depending on the parameters. 
The steps following the first one is due to the development of a vibrational instability and a 
transition into the associated charge transfer regime. The obtained behavior of the I-V-curves 
is in reasonably good agreement with the experimental data 0. Best fit to the published 
experimental I-V-curves is obtained for an asymmetric coupling to the leads {Th/Tl ~ 9). 
When the ratio Th/Tl decreases both the vibration-induced current jumps and a high- voltage 
slope of the I-V-curve increase deviating from the experimental data. 

An alternative theoretical description of the experiment based on a photon-assisted tunneling 
like picture Il3] also shows reasonable agreement with the experimental data M . Therefore it 
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Fig. 3 - I-V curves for different values of the parameter x which characterizes the strength of the 
electric field between the leads for a given voltage: wo = 5meV, T = 0.13 meV, eo = 6meV, 
and r = 2.3 //eV. Best fit to Q is obtained for an asymmetric coupling to the leads; here we use 
Vr/Vl = 9. 



is difficult to conclude with certainty that it is the above instability which is responsible for 
the specific features observed in the experiment. Further experiments are needed to clarify 
the picture. Such experiments may involve an investigation of the charge fluctuation on the 
grain or of current noise. 

Conclusion. - We have studied the effect of a couphng between coherent tunneling of 
electrons through a single quantized energy level in the central island — or dot — of a single 
electron transistor and vibrations of the dot. We have found that for bias voltages exceeding 
a certain critical value a dynamical instability occurs and mechanical vibrations of the center 
of mass of the dot develop into a stable limit cycle. The effect of this vibrations on the current 
through the system were also studied. I-V characteristics calculated in our model were found 
to be in a reasonably good agreement with recent experimental results of Park et al. Q for 
the single Ceo-molecule transistor. 
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